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Introduction 


EXCEDE  m  (Excitation  by  Electron  Deposition)  is  an  atmospheric  energy  deposition  rocket 
experiment  which  was  launched  on  27  April  1990  from  White  Sands  Missile  Range.  Included 
as  part  of  the  EXCEDE  II  sensor  module  ("daughter")  payload  were  several  photometers  with 
a  mechanically  scanned  mirror.  The  purpose  of  these  instruments  was  to  obtain  spatial  profiles 
of  the  beam  intensity  at  a  variety  of  wavelengths  and  at  a  fixed  angular  distance  along  the  beam 
from  the  accelerator  module.  One  of  these  photometers  was  filtered  to  respond  to  light  at  3914 
A. 


Emissions  at  3914  A  result  from  the  excitation  of  the  (0,0)  first  negative  band  of  N*2 
[denoted  N2+  (IN)  Bom  and  Zipf,  1970].  The  excited  state  produced  by  electron  impact  is  B2 
which  has  a  lifetime  of  70  ns  [Jones,  1971]  and  decays  to  X3  E+g.  This  excitation  has  a 
threshold  energy  of  18.8  eV,  a  peak  cross  section  at  100  eV,  and  a  slow  decrease  in  the  cross 
section  toward  higher  energies  [Bom  and  Zipf,  1970].  These  factors  make  emissions  at  this 
wavelength  an  excellent  indicator  of  the  spatial  extent  of  energy  deposition  by  the  primary  beam. 
Monitoring  of  this  line  has  been  a  standard  diagnostic  tool  in  auroral  observations. 

Experiment  Goals 

The  primary  goal  for  the  scanning  photometer  investigations  was  to  obtain  the  spatial 
distribution  of  energy  deposited  in  the  upper  atmosphere  by  the  EXCEDE  m  beam.  The 
technique  employed  the  3914  A  emission  as  a  deposition  diagnostic  (also  used  in  the 
characterization  of  natural  aurora).  The  narrow-band  photometer  (  — 13  A  FWHM)  was  scanned 
from  -15°  (ahead  of  the  beam)  to  +30°  (trailing  the  beam)  across  the  nominal  beam  location 
at  —  1  scan  per  second.  The  measured  field  ofview  for  this  photometer  is  —0.1°.  The  viewing 
geometry  is  shown  schematically  in  Figure  1.  A  variety  of  other  infrared,  visible,  and 
ultraviolet  photometers  and  spectrometers  were  also  used  to  measure  intensities  at  the  beam  and 
in  the  afterglow  region.  The  results  of  these  studies  will  be  presented  elsewhere. 

The  scanning  photometer  experiment  package  was  designed  and  assembled  at  Visidyne  and 
consists  of  four  photometers  with  narrow  band  filters,  a  motor-driven  scanning  mirror,  the  optics 
assembly  and  associated  electronics.  A  schematic  drawing  of  the  instrument  is  shown  in  Figure 
2. 


The  3914  A  photometer  consisted  of  a  EMI  54  IN  photomultiplier  tube  (PMT)  in  conjunction 
with  an  aperture,  a  lens,  and  a  narrow  band  filter.  The  sensor  was  limited  to  maximum  counting 
rates  of  —2x10s  counts/second  and  had  dark  noise  counts  of  38  per  second.  As  mentioned 
above,  the  narrow  band  filters  have  characteristic  widths  (FWHM)  of  — 13  A.  Recent 
measurements  of  these  same  transmission  curves  show  that  the  spectral  response  of  the  filter  has 
remained  stable  two  years  after  the  flight. 

In  spite  of  the  relatively  high  emission  intensity  of  the  beam,  the  number  of  counts  in  each 
measurement  bin  are  relatively  low  during  many  of  the  scans  near  apogee  because  the  integration 
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Figure  1.:  Viewing  geometry  of  the  sensor  module  with  respect  to  the  accelerator  module. 
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SCANNING  MIRROR 


Figure  Z:  Schematic  drawing  of  the  EXCEDE  m  scanning  photometer. 
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time  per  measurement  is  short  (1.33  ms).  The  instrument  count  rate  was  set  by  the  need  for 
minimizing  the  uncertainties  due  to  counting  statistics  while  maximizing  the  spatial  resolution 
across  the  beam. 

In  addition  to  the  photometer  centered  on  3914  A,  the  unit  also  included  narrow-band 
photometers  centered  on  3805  A  (Ni(2P)  p'=0,v"~2),  5577  A  (*S0  -  'D2  "auroral  green  line"), 
and  2761  A  (N2  Vegard  Kaplan  (0,6)  transition).  Both  of  the  excited  states  for  the  latter  two 
emissions  are  long-lived.  The  emissions  measured  by  these  photometers  trace  out  the  evolution 
of  the  region  into  which  the  beam  energy  is  deposited  and  will  be  discussed  elsewhere. 

Experimental  Results 

t 

There  were  170  scans  made  by  the  3914  A  photometer  between  133  s  (after  launch)  and  300 
seconds,  including  90  scans  with  the  electron  beam  on.  We  have  analyzed  the  beam-on  scans 
in  terms  of  a  generalized  Gaussian  function  (described  below)  plus  a  background  count  rate. 
The  beam-off  scans  were  used  only  to  determine  a  background  photon  count  rate. 

In  the  3914  A  data  there  is  an  asymmetry  in  the  scan  profiles  early  in  the  flight  which 
develops  into  a  second  completely  separated  peak  by  the  time  of  165  seconds.  The  secondary 
peak  always  appears  on  the  side  of  the  primary  peak  that  is  "downstream"  with  respect  to  the 
motion  of  the  two  modules  through  the  ionosphere.  The  ratios  of  the  areas  under  the  two 
intensity  curves  varies  from  -  5  to  - 10% ;  the  widths  are  similar  (although  less  well-determined 
for  the  secondary  peak).  The  angular  separation  of  the  two  remains  relatively  constant  at  1.75°. 
This  latter  constancy  of  angle  suggested  that  the  second  peak  could  be  an  instrumental  artifact 
(F.  Bien,  private  communication). 

Subsequent  to  this  suggestion,  the  flight  photometer  unit  was  set  up  in  the  laboratory  to  test 
this  hypothesis  (the  instrument  payload  on  the  Sensor  Module  was  recovered  after  successful 
parachute  deployment).  A  calibrated  NIST  source  lamp  employing  a  halogen-quartz  envelope 
and  coiled-coil  filament  was  used  to  produce  an  extended  source  by  illuminating  a  diffusion 
screen.  The  screen  was  then  masked  to  look  like  a  spatial  line  source  similar  to  the  EXCEDE 
ID  beam  (full  width  —0.75°).  The  source  was  then  scanned  with  the  flight  unit  mirror  with  the 
intensity  from  the  3914  A  photometer  recorded.  A  schematic  of  the  laboratory  setup  is  shown 
in  Figure  3.  This  laboratory  simulation  revealed  a  secondary  peak  with  an  integrated  area  of 
10%  of  the  primary  peak  and  at  an  angular  separation  of  1.68°.  Further  investigation  has  shown 
that  the  extra  peak  in  the  3914  A  data  can  be  traced  to  a  reflection  due  to  the  narrow-band  filter 
assembly. 

The  effects  of  the  reflection  on  the  flight  data  are  minimal  and  do  not  compromise  any  of 
the  results.  Because  the  reflection  always  occurs  on  the  same  side  of  the  scan,  the  leading  edge 
of  all  the  scans  are  not  affected.  The  data  are  fit  to  a  symmetric  function  and  as  can  be  seen 
in  Figure  4  the  unaffected  edge  of  the  fit  is  in  excellent  agreement  with  the  data.  Where 
possible,  we  have  fit  the  data  to  a  two-peak  intensity  distribution  and  then  analyzed  the 
parameters  of  the  primary  peak  alone  to  determine  the  energy  deposition  properties  of  the  beam. 
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LABORATORY  SETUP 


Figure  3.i  Schematic  of  post  calibration  laboratory  setup. 


5 


Unfolding  the  Deposition  Profile.  Let  p(r)  be  the  volume  emissivity  of  the  beam,  where  r 
is  the  distance  from  the  beam  axis  to  the  differential  volume  of  emission  in  cylindrical 
coordinates.  Here  we  assume  that  the  emission  pattern  (with  the  secondary  peak  eliminated)  has 
cylindrical  symmetry.  Let  the  x-axis  be  defined  by  a  line  along  the  instantaneous  view  direction 
of  the  photometer,  let  the  y-axis  be  along  the  direction  perpendicular  to  the  x-axis  and  through 
the  nominal  center  of  the  beam,  and  let  the  z-axis  be  coplanar  with  the  symmetry  axis  of  the 
beam.  We  assume  that  the  observer  is  located  a  distance  r  from  the  nominal  beam  center  and 
that  the  x-axis  makes  an  angle  6  to  the  beam  axis;  the  geometry  is  indicated  in  Figure  5.  The 
distance  y  is  then  related  to  the  photometer  mirror  scan  angle  4>  by 

y  -  R  sin  4>  ^ 

If  E'(y)  is  the  measured  intensity,  and  E  is  the  intensity  that  would  be  observed  perpendicular 
to  the  beam  axis,  then  we  have 

E(y)  =  E'(y)  sin  6 


and 


(3) 


where  i2  =  x2  +  y2.  Equivalently 


£(y) 


7  dr 

i  yjj 2  -  y 2 


(4) 


where  we  have  made  use  of  the  assumed  cylindrical  symmetry  by  writing  p(r)  for  p(x,y). 


The  volume  emissivity  is  given  by  the  solution  of  eqn.  (4),  which  is  a  Volterra  integral 
equation  first  studied  by  Abel  [Courant  and  Hilbert ,  1953].  The  solution  leads  to  the  Abel 
transform,  which  is  encountered  in  unfolding  projected  luminosity  functions  in  a  variety  of 
settings  [BraceweU,  1956,  Press  et  al.,  1976].  The  emissivity  p  is  given  by  the  inverse  Abel 
transform  [e.g.,  Bracewell ,  1956] 


QS 


is 


(5) 


Generalized  Gaussian  Distribution.  The  statistical  noise  in  the  data  precludes  the  direct 
application  of  the  inversion  in  Eqn.  (5).  What  can  be  inverted  is  an  analytic  form  for  E(y) 
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which  includes  a  sufficient  number  of  parameters  to  characterize  the  data  well.  These 
parameters  (plus  their  statistical  uncertainties)  can  then  be  obtained  by  solving  the  "forward 
problem"  in  a  least  squares  sense  (e.g.,  [Bevington,  1969]). 

As  a  first  cut  to  the  interpretation  of  the  3914  A  data  we  have  found  a  convenient  analytic 
function  for  fitting  the  intensity  scans  to  be  a  "generalized  gaussian"  defined  by  [Tarantola, 
1987] 


2,09 


Az  sin0 


pi-vp 

— - - exp 

2<r  T(l/p)  F 
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with 

♦0 

|  2,09  dy  =  Az  sin0 


(7) 


and  where  y0  is  the  center  of  the  distribution,  a  is  a  measure  of  die  width,  and  p  is  the  power 
in  the  generalized  Gaussian  ( for  p=2  the  distribution  reduces  to  a  simple  Gaussian).  By  using 
a  generalized  Gaussian  for  Z(y)  measured  from  the  symmetry  point,  i.e.,  putting  y„  =  0  we 
obtain 


P(r) 


1  Az  sin0 

x  Ta/pT 
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(pUp  a)2 
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where 
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p1,p  cr 


and 
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or,  equivalently, 


fp  (a)  =  p  a1  •  Up  |  e-aa**'  CQstf~lv  dv 


(ID 


9 


In  general  £  (a)  cannot  be  evaluated  in  terms  of  simpler  functions;  however,  there  are  simple 
forms  for  the  following  special  cases 

A  («)  *  K0(a) 

fi  («)  ■  e" 
ft  (0)  -r(l  -  Up)  for  p  >  1  (14) 

3914  A  Individual  Scans.  The  90  scans  made  by  the  3914  A  photometer  with  the  electron 
beam  on  have  been  fit  to  a  generalized  gaussian  with  the  parameters:  AE  *  integrated  intensity 
(Megarayleigh  -  meter);  y0  *  colter  location  (meter);  a  *  width  (meter);  p  ■  generalized 
exponent;  and  I0  =  constant  background.  The  beam-off  scans  were  averaged  to  compute  a  noise 
level.  Plots  of  the  five  fitted  parameters  are  displayed  in  Figures  6  -  10  as  a  functions  of  time. 

Figure  6  is  a  plot  of  the  amplitude,  As,  in  units  of  MR-meters.  Ac  is  equal  to  the  area 
under  the  generalized  gaussian  and  represents  the  total  number  of  3914  A  photons  emitted  from 
the  scanned  slice  of  the  beam.  The  distribution  of  amplitudes  is  approximately  symmetric  about 
195  seconds  (apogee)  where  the  values  are  close  to  minimum.  The  3914  A  intensity  is 
characterized  by  the  atmospheric  density  being  largest  at  times  corresponding  to  lower  altitudes 
and  smallest  at  higher  altitudes.  The  two  clusters  of  low  points  (centered  at  208  and  215 
seconds)  correspond  to  gun  cycles  13  and  14  where  the  gun  currents  were  reduced  due  to  load 
faults  in  the  accelerator  system. 

The  generalized  exponent,  p,  is  plotted  in  Figure  7.  From  the  plot,  we  can  see  that  a  typical 
value  for  p  lies  between  1.3  and  1.8.  The  value  of  these  exponents  infers  that  the  beam  is  more 
peaked  at  its  center  than  would  be  the  case  for  a  standard  gaussian  distribution  (i.e.  p=2). 

The  location  of  the  center  of  the  beam  is  plotted  in  Figure  8  in  units  of  degrees  and  meters. 
A  value  of  zero  indicates  that  the  center  of  the  interferometer  FOV  (from  which  all  instrument 
directions  are  referenced)  is  centered  on  the  beam.  Positive  (negative)  values  of  y„  correspond 
to  the  center  of  the  interferometer  being  aimed  ahead  of  (behind)  of  the  beam.  The  derived 
centroids  indicate  that  the  interferometer  is  initially  aimed  ahead  of  the  beam  by  a  few  degrees 
and  becomes  more  nearly  centered  up  to  -250  seconds. 

The  width  of  the  generalized  Gaussian,  <r  (see  eqn.  6),  is  also  plotted  in  units  of  both  degrees 
and  meters  in  Figure  9.  The  value  of  a  stays  relatively  constant  at  -  2.3  meters  over  the  entire 
upleg  of  the  flight  and  then  steadily  increases  to  over  4  meters  at  250  seconds.  Note,  however, 
that  in  terms  of  angular  size  a  remains  essentially  constant  at  0.35  degrees  through  the  downleg. 
However,  because  the  generalized  exponent  can  also  change,  the  width  of  the  beam  is  better 
characterized  by  the  RMS  diameter  described  in  the  following  section. 
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Photometer  Fit  Parameters 
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3914  Scanning  Photometer  Fit  Parameters  -  Generalized  Exponent 


3914  Scanning  Photometer  Fit  Parameters 
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parameters  -  centroid  (In  units  of  meters  and  degrees) 


3914  Scanning  Photometer  Fit  Parameters 
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Mission  Elapsed  Time  (Seconds) 

3914  A  scanning  photometer  fit  parameters  -  generalized  width  (In  units  of  meters  and  degrees). 


3914  Scanning  Photometer  Fit  Parameters  -  Background 
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Mission  Elapsed  Time  (seconds) 

I  A  scanning  photometer  fit  parameters  -  background, 


The  fitted  background  intensity  is  plotted  in  Figure  10  and  die  background  with  die  beam  on 
tends  to  be  proportional  to  the  integrated  3914  A  signal  characterized  by  Ar  throughout  the 
flight  This  is  to  be  compared  with  the  beam  off  scans  which  have  a  nearly  constant  intensity 
(—0.004  MR)  independent  of  time  or  altitude.  The  background  at  3914  A  is  -0.006  MR  near 
apogee  which  corresponds  to  -196  counts/seconds  or  -S  times  the  dark  noise  of  die 
photomultiplier  detector.  This  difference  may  be  caused  by  scattered  light  inside  the  photometer 
or  may  be  an  artifact  of  the  fit. 


Transmission.  The  spectral  response  of  the  narrow  band  filter  used  in  the  photometer  is 
gaussian  in  shape  whereas  the  spectral  shape  of  the  3914  A  emission  is  asymmetric  with  a  short 
wavelength  tail  resulting  from  its  band  structure.  Consequently,  the  filter  does  not  cover  the 
entire  spectral  region  of  the  3914  A  emission  and  the  spectral  "efficiency"  must  be  calculated. 
Figure  11  shows  a  plot  of  the  transmission  of  the  filter  and  the  spectral  scan  at  205  seconds 
(114.5  km).  The  spectral  efficiency  e  is  calculated  by 


(15) 


where  z  is  the  altitude  and  IMU(X)  is  the  spectral  band  profile.  If  had  the  shape  of  a  top 

hat  and  encompassed  the  entire  3914  peak,  e  would  be  equal  to  1. 


Fourteen  visible  spectrometer  scans  were  used  to  determine  the  value  of  c.  The  integral  over 
the  N2*  band  was  carried  out  numerically  by  amply  summing  the  spectral  bins  between  die 
wavelength  region  from  3870  A  to  3940  A.  The  value  for  e  showed  no  significant  systematic 
variations  over  die  14  independent  evaluation  and  were  therefore  averaged  together  to  obtain  a 
value  of  0.546  ±  0.012.  This  represents  the  fraction  of  the  Nj+  emission  that  was  observed; 
die  entire  3914  A  data  was  then  scaled  by  1/e. 

Comparison  with  Photographic  Data.  Three  data  scans  from  the  3914  A  scanning 
photometer  were  compared  with  similar  scans  derived  from  a  photographic  analysis  made  by 
Photometries,  Inc..  Scans  intended  to  represent  the  overall  flight  were  chosen  at  the  following 
times:  152  seconds  (pulse  5),  201.7  seconds  (pulse  12),  and  258.5  seconds  (pulse  20).  The 
results  are  shown  in  Figures  12, 13,  and  14.  The  solid  line  is  the  fit  to  the  scanning  photometer 
data  by  a  generalized  gaussian  function;  the  points  are  data  points  from  a  micro-densitometer 
scan  of  the  photographic  data.  The  amplitudes  are  relative  and  were  scaled  for  the  comparison. 
The  agreement  of  the  shapes  at  152  and  201.7  seconds  is  impressive  and  provides  confidence 
in  the  measurement  of  the  energy  deposition.  The  disagreement  in  the  two  data  sets  at  258.5 
seconds  is  likely  due  to  the  two  instruments  viewing  different  locations  on  a  rapidly  changing 
beam.  It  is  also  noted  that  the  amplitudes  from  the  two  measurements  have  been  reconciled  and 
agree  to  better  than  15  percent.  Additional  analysis  of  the  photographic  data  is  being  done  and 
will  be  reported  elsewhere. 
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Response  of  3914  A  Filter 


Figure  11. 


EXCEDE  III  -  3914  Scanning  Photometer  :  152  Seconds 
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Figure  12.  Comparison  of  scanning  photometer  data  with  photographic  data  at  152  s. 


EXCEDE  III  -  3914  Scanning  Photometer  :  201.7  Seconds 
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Figure  19.  Comparison  of  scanning  photometer  data  with  photographic  data  at  201.7  s. 


EXCEDE  III  -  3914  Scanning  Photometer  :  258.5  Seconds 
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Figure  14.  Comparison  of  scanning  photometer  data  with  photographic  data  at  258.5  s. 


Analysis 


We  have  obtained  good  fits  for  -  40  scans  prior  to  apogee  and  for  -  40  past  apogee  (apogee 
was  115  km  and  occurred  195  s  into  the  flight). 

One  possible  way  to  characterize  the  measured  width  of  the  beam  is  to  compute  the  RMS 
radius.  This  characterization  is  independent  of  the  exponent  (  other  characterizations  such  as 
FWHM  are  not)  and  represents  the  point  where  half  of  the  integrated  power  is  within  that 
distance  and  half  of  the  power  lies  outside  of  that  dimension. 

The  RMS  radius  is  defined  as  the  second  moment  of  the  distribution  and  is  described  by  the 
expression  (see  Appendix  A.) 

o§ 

[  p(r)  r3  dr 

<r2>  m  2 _  =  0*  Qp2Jf)  tt® 

-  ni/p) 

f  p(r)  r  dr 


The  RMS  radius  is  plotted  as  a  function  of  altitude  in  Figure  15  and  remained  relatively 
constant  above  ~  102  km,  ranging  between  —  4  and  5  meters.  The  increase  in  the  beam  radius 
seen  at  lower  altitudes  is  due  to  viewing  the  beam  at  a  larger  fraction  of  the  practical  range 
during  this  portion  of  the  flight.  The  beam  diameter  is  increased  because  of  the  larger  number 
of  collisions. 

The  atmospheric  dosing  from  the  EXCEDE  m  beam  has  been  calculated.  The  beam  is  an 
extended  source  and  consequently  the  dose  is  different  at  different  points  across  the  beam.  The 
dose  is  related  to  the  energy  deposition  by 

D(y,t )  * 

where  E(y)  is  the  energy  deposition  as  determined  from  the  3914  A  scanning  photometer  and 
vx  is  the  velocity  component  of  the  accelerator  module  transverse  to  the  local  magnetic  field. 

The  number  of  detected  3914  A  photons  can  be  related  to  the  energy  deposition  by  noting 
that  for  every  3914  A  photon  there  are  — 14. 1  N2+  ions  created  [Borst  and  Zipf,  1970] 
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EXCEDE  III  —  RMS  Beam  Radius 
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Figure  15.  Plot  of  RMS  radius  derived  from  3914  A  scanning  photometer  data 


(assuming  the  electron  energies  are  in  excess  of  -SO  eV).  The  fraction  of  ion  pairs  that 
produce  N2+  ions  is  given  by  [Rees  and  Jones,  1973] 


ionization  fraction  - 


0.92  [A£] 

1.15  [A£]+ 1.5  [02*1+0.56  [M?*] 


and  each  ion  pair  has  associated  with  it  the  local  deposition  of  —  35  eV  of  energy. 

An  average  dosing  rate  can  be  calculated  assuming  the  RMS  beam  diameter  as  a 
characteristic  width  for  the  beam.  The  average  dosing  rate  is  given  by 

<b»>,  -  ±  i*Mdy  m 

’  <y  > 

This  quantity  is  easily  calculated  by  noting  that  the  integral  over  the  distribution  E(y),  is  equal 
to  the  fit  parameter  Ac. 


At  sin  6  = 


=  |  Z'p  (y)  dy  =  2r  |  p  (r)  r  dr 


The  average  dose  is  calculated  by  dividing  AE  by  vx  and  the  RMS  diameter.  A  plot  of  the 
average  and  peak  dosing  rates  is  shown  in  Figure  16. 

The  maximum  dose  is  reached  on  upleg  at  - 102  km  and  is  in  excess  of  1013  eV/cm3.  The 
value  of  the  dose  gets  smaller  all  the  way  up  to  apogee  where  it  is  a  minimum  (-3.5  x  1010 
eV/cm3)  and  then  increases  as  the  altitude  decreases. 

Caution  is  required  when  considering  the  maximum  dose  values.  Let  (dP/dx),  be  the  energy 
deposited  per  second  along  the  beam  at  a  distance  s  from  the  accelerator.  The  assumption  is 
made  for  the  above  calculations  that  the  value  for  dP/dx  is  approximately  constant  close  to  the 
scanned  region.  This  is  a  very  good  approximation  for  the  majority  of  the  flight  where  the  dose 
time  is  short  and  the  accelerator  module  moves  very  little  in  the  beam  direction.  Under  these 
conditions  the  scanned  region  is  dosed  by  the  portion  of  the  beam  approximately  a  distance  s0 
along  the  beam.  This  is  not  true  where  the  dose  time  is  a  maximum.  In  that  case  almost  all  of 
the  motion  is  along  the  beam  direction  and  therefore  the  irradiated  volume  is  exposed  to  an 
energy  deposition  rate,  (dP/dx),  at  many  different  values  of  s  on  the  beam  which  varies 
significantly. 
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EXCEDE  III  -  Average  and  Peak  Dose 
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average  and  peak  dose  derived  from  3914  A  scanning  photometer  data 


Similarly,  an  average  surface  brightness  for  the  3914  A  emission  can  be  calculated  by  the 
following  expression 


sing  he 
4t  2 rjats  A 


(21) 


Catalog.  A  table  of  all  fit  parameters  have  been  compiled  for  all  scans  and  is  presented  in 
Appendix  B.  in  addition  to  the  relevant  fit  parameters  and  the  associated  errors,  the  RMS 
diameter  and  dosing  for  each  scan  has  been  included.  Note  that  although  a  secondary  gaussian 
(p=2)  was  fit  where  possible,  the  primary  peak  is  taken  to  describe  the  beam  deposition. 

Plots  of  the  data  and  the  corresponding  fit  are  shown  in  Appendix  C.  Each  scan  is  displayed 
on  a  linear  scale  in  the  upper  plot  and  a  logarithmic  scale  in  the  lower  plot.  Plots  with  unusually 
large  values  for  p  (e.g.  scan  number  1.)  correspond  to  incomplete  scans  where  die  beam  was 
cither  turned  on  or  off  during  the  data  collection. 

A  floppy  disk  containing  the  Appendix  B  and  the  data  set  presented  in  Appendix  C  has  also 
been  produced. 

Calculating  the  Abel  Transform.  In  many  cases  the  integrated  line  scan  across  the  beam,  A E 
is  adequate  for  a  particular  application  (e.g.  calculating  an  auroral  efficiency).  For  others, 
however,  the  derived  radial  distribution  is  necessary.  In  the  latter  case  the  corresponding 
inverse  Abel  transform  must  be  calculated.  A  method  is  described  below  for  extracting  the 
value  of  the  radial  distribution  at  a  point  r  for  a  given  scan.  The  following  approach  reduces 
the  expression  to  a  dimensionless  quantity  described  by  a  family  of  curves. 

Each  spatial  scan  from  the  3914  A  scanning  photometer  has  been  fit  to  a  generalized  gaussian 
distribution  and  can  be  inverted  using  the  inverse  Abel  transform  to  obtain  the  radial  distribution. 
The  radial  distribution  can  be  expressed  in  the  generic  form 

p(f)  *  sind  P(-,p)  (22) 

2  xa2  o 


The  parameters  Ac,  a,  p,  sind,  and  the  generalized  exponent  p  are  listed  in  Appendix  B  for  each 
scan  along  with  other  parameters. 
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As  an  example  of  this  method  consider  pulse  4  att=  136.83  seconds.  What  is  the  value  of 
p  at  a  distance  of  5  meters  from  the  beam  center?  The  corresponding  value  of  the  fit 
parameters  at  t=136  seconds  are 


Ac  *  198.86  MR-meter 
a  =  2.11  meters 
p  *  1.40 
sin0  »  0.69 


r/a  *  2.37 


Figure  17  shows  a  family  of  curves  of  P(r/<r,p)  for  different  values  of  p;  specific  values  of 
P(r/«r,p)  can  be  obtained  from  this  plot.  The  value  of  P(2.37, 1.40)  is  found  to  be  equal  to  0.073 
from  Figure  17.  The  value  for  e  is 


p<2.37) 


(198.86) 

(2r)  (2.11  x  102  cm)2 


(0.69)  (0.073)  MR -meters  x  100  cm/meter 


(23) 


■  3.58  x  10*  photons! cm1 -s 

Note  that  the  units  are  Megarayeigh-meters/cm2  and  have  been  converted  to  photons/cm3- s. 
The  user  may  ultimately  find  it  easier  to  calculate  P(r/a,p)  using  the  formulation 

P(L,p)  =  al~Up  |  "  cosh''1?  dv  (24) 


where 


a  -  1  (Ly  (25) 

P  <* 

which  is  computationally  better  behaved  at  r=0. 

Figures  18  and  19  are  plots  of  the  energy  deposition  from  the  beam  for  the  above  fit 
parameters  and  are  intended  to  characterize  the  overall  beam  profile.  The  Abel  Transform  has 
been  calculated  for  the  values  p=1.40  and  <r=2.11  meters  (see  scan  number  46)  and  the 
resulting  spatial  distribution  has  been  plotted  as  an  image  in  Figure  18.  Figure  19  is  a  2- 
dimensional  picture  of  the  same  cross  sectional  distribution. 

Summary 

The  energy  deposition  from  the  EXCEDE  m  beam  was  monitored  by  observing  the  N^flN) 
emission  at  3914  A  with  a  narrowband  filter  photometer  that  scanned  across  the  beam  each 
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Inverse  Abel  Transform  of  Generalized  Gaussian  Distribution 
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Figure  17.  Plot  of  the  family  of 


Typical  EXCEDE  III  Beam  Profile 


scanning  photometer  measurement. 


second.  Ninety  (90)  scans  were  collected  with  the  beam  on  and  were  used  to  characterize  the 
spatial  deposition  of  the  beam.  Quality  fits  to  the  data  were  obtained  using  generalized  Gaussian 
distributions.  These  distributions  were  then  inverted  by  calculating  the  inverse  Abel  transform 
to  obtain  the  corresponding  radial  distributions  of  the  energy  deposition. 

The  RMS  radius  for  each  scan  of  the  beam  has  been  calculated  along  with  the  average  and 
peak  dose  which  is  in  excess  of  10u  eV-cm'3  near  maximum  dose. 
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Appendix  A:  Moments  of  the  Distributions 


It  is  convenient  to  describe  the  deposition  profile  in  terms  of  the  spatial  moments  of  the 
profile.  By  making  use  of  the  definition  of  the  Abel  transform  and  changing  the  order  of 
integration,  several  useful  relationships  can  be  derived.  For 


M. 


(26) 


we  obtain 


M. 


r(iLLi) 

— £—  (<r  PVryt 

T(±) 

P 


(27) 


so  that,  in  particular.  Mo  =  1/2  and  Mp  =  1/2  a9.  By  making  use  of  Eqn.  (4)  and  using  that 
with  M,  as  defined  by  Eqn.  (25),  we  can  obtain 


(££!) 

TV*  *  2) 

2 


f  p  (r)  r"  * 1  dr 


(28) 


% 

It  is  now  trivial  to  show  that  the  rms  radius  of  the  deposition  (i2  averaged  over  the  emissivity) 
is  then  given  by 


<r*> 


p(r)  r3  dr 

-  =  a2  (2 p2") 

p(r)  r  dr 


TQlp) 

T(Up) 


(29) 
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